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SUMMARY 


A rapid computer program for the calculation of the design and the off-design 
performance of a turbojet or turbofan engine is described. The program is 
designed to be run in a stand-alone mode or operated within a larger program 
such as ACSYNT, an ^ir£raft synth esis program for preliminary design, which is 
in use at Ames Research Center. The computation is based on a simplified one- 
dimensional gas turbine cycle. Each component in the engine is modeled thermo 
dynamically. The component efficiencies used in the thermodynamic modeling 
are scaled for the off -design conditions from input design point values using 
empirical trends which are included in the computer code. The engine cycle 
program is capable of producing reasonable engine performance prediction with 
a minimum of computer execute time. The current computer execute (CPU) time 
on the IBM 360/67 for one Mach number, one altitude, and one power setting is 
about 0.1 seconds. 

The principal assumption used in the calculation is that the compressor is 
operated along a line of maximum adiabatic efficiency on the compressor map. 
The fluid properties are computed for the combustion mixture, but dissocation 
is not included. The procedure included in the program is only for the com- 
bustion of JP-4, methane, or hydrogen; however, modifications are possible 
which would allow the calculation for any hydrogen, carbon, or hydrocarbon. 
Sample input and output listings are included. 


INTRODUCTION 


The design of advanced aircraft requires the consideration of many parameter 
variations in order to arrive at an optimum design for a particular mission or 
group of missions. These design studies are traditionally done using many 
large, separate computing programs that consider the various disciplines in a 
noninteractive sense. However, recent rapid advances in computer technology 
allow use of the computer to automate the early conceptual design of airborne 
systems. The desire to use the computer for aircraft conceptual design 
generates the need for a rapid, accurate means of calculating the performance 
of the propulsion system. This report is oriented toward satisfying this need 



The most detailed qas turbine cvcle analyses for desiqn and off-desiqn 
performance are 6ENENG ( GENe ralized ENGi ne) {refs. 1 and 2) which was developed 
from SMOTE (Simulation of X^^i^bofan Engine) (ref. 3) and NNEP (IMavy-NASA £ngine 
Program) which was derived from NEPCOMP (refs. 4 and 5, respectiveliO- These 
programs perform a detailed flow and performance match of the compressor and 
turbine for each desired operating point. This matching requires the use of 
detailed compressor and turbine performance maps and a means of matching the 
compressor and turbine performance at each operating point. These programs 
run in slightly less than 10 seconds per Mach number, altitude and power setting 
on an IBM 360/67. 

Waters (ref. 6) has published a description of a design or off-design gas 
turbine cycle analysis. This analysis has an empirically determined effect 
of compressor speed variation upon off-design compressor pressure ratio, but it 
does not consider the effect of off-design operation on the various component 
efficiencies. 

Vanco (ref. 7) has published a design-point cycle calculation. This work does 
not consider the off-design case, nor does it have the built-in capability to 
size the engines for a particular mission requirement. 

The computer code presented in this report can compute the design or off-design 
performance of a one or two spool turbojet or turbofan engine, with or without 
afterburning or ductburning. It uses a compressor model which empirically 
estimates the off-design pressure ratio and adiabatic efficiency of the com- 
pressor. The efficiencies of the other engine components are also modeled 
empirically. The gas properties are approximated using the polynomial approxi- 
mations suggested by Vanco (ref. 7). The present program uses a model for the 
off-design fan pressure ratio which eliminates the cycle iterations required 
for matching the engine. Use of the presented program increase the computa- 
tional speed by at least a factor of 10 when compared with other programs which 
use cycle iterations. This computational speed is achieved at the expense of 
decreased accuracy. 

A comparison of the results of the computer program with the performance of 
existing engines is included. 


SYMBOLS 

A flow area (meters^, feet^) 

BR bypass ratio 

Cp specific heat at constant pressure (Joules/kg°k, BTU/lbm°R) 

dp derivative of total pressure (Newton/meter^, Ibf/ft^) 

f fuel air ratio by mass 

g gravitational constant (meters/sec, ft/sec^) 
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h 

HVF 

J 

M 

N 

P 

PRCD 

PRFD 

R 

RPM 

s 

SFC 

T 

V 

W 

AS 

A({) 

n 

SUBSCRIPTS: 

1 

2 

3 

4 

5 

6 
7 


specific total enthalphy (Joules/kg, BTU/lbm) 

fuel lower heating value (Joules/kg, BTU/lbm) 

conversion factor (1.0 for SI units, 778 ft Ibf/BTU) 

molecular weight 

rotational speed (RPM) 

total pressure (Newtons/m , Ibf/ft ) 

compressor design pressure ratio 

fan design pressure ratio 

universal gas constant (Joule/kgmol°K, ft lbf/lbmo°R) 

revolutions per minute 

entropy 

specific fuel consumption 

total temperature (°K, °R) 

velocity (meters/sec, ft/sec) 

mass flow (kg/sec, Ibm/sec) 

incremental enthalpy (Joules/kg, BTU/lbm) 

incremental entropy 

incremental entropy function 

efficiency 


station numbers (see figure 1) 
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a 

initial 

A/B 

afterburner 

air 

air 

ambient 

ambient 

b 

final 

B 

burner 

comp 

compressor 

design 

design condition 

D/B 

duct burner 

f 

fuel 

i 

isentropic (or ideal) 

in 

into considered component 

nozzle 

nozzle 

0 

free-stream 

out 

out of considered component 

ov 

overal 1 

prim 

primary 

product 

products of combustion 

pro 

propulsive 

R 

reference 

s 

static 

sec 

secondary (duct) 

th 

thermal 

t 

turbine 


CYCLE DESCRIPTION 


A schematic diagram of a genera! gas turbine engine cycle is shown in Figure 1. 
The air enters the inlet with a resultant loss in stagnation pressure which is 
typically a function of flight Mach number. The airflow then enters the fan 
and/or outer compressor where it may be divided into the primary flow (the flow 
through the primary gas turbine) and, if the engine is a turbofan, the second- 
ary or bypass flow. The primary airflow is then routed through the inner com- 
pressor for additional compression. A small amount of the primary flow is 
directed around the combustor to provide cooling for the turbines. The 
remainder of the primary flow enters the combustor where fuel is added and 
combustion occurs. The combustors are designed for fuel -lean conditions (less 
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fuel than required for stoichemetric burning) in order that the gas temperature 
entering the turbine does not exceed the structural or material limitations of 
the turbine. Many studies have been made on methods of cooling these turbines 
and quite involved techniques are used in modern gas turbine engines. Current 
cooling technology allows the turbine inlet temperature to go as high as 
1389°K to 1778°K (2500°R to 3200°R). 

The primary airflow, including the turbine cooling air that was routed around 
the combustor, enters the inner turbine which is connected to the inner com- 
pressor with a shaft along the centerline of the engine. This shaft transfers 
the energy required to drive the inner compressor from the inner turbine to the 
inner compressor. The combination of compressor, shaft, and turbine is usually 
called a spool. The primary gas flow then enters the outer turbine which is 
connected via a shaft running through the hollow inner turbine shaft to the 
outer compressor-fan combination. This shaft supplies the energy required to 
drive the outer compressor-fan combination. The combination of outer compressor- 
fan, outer turbine and connecting shaft is called the outer spool. The fan 
may or may not have the same number of stages (rows of blades) as the outer 
compressor. In a three spool engine, the fan is driven by a separate turbine. 
Many engines do not have a fan but have an inner and outer compressor (a two 
spool turbojet). The gas flow leaving the outer turbine enters the afterburner 
section where additional fuel can be added and burned. The afterburner may or 
may not be present, depending on the operational requirements of the engine. In 
general, the afterburner supplies a large amount of additional thrust (about 1.5 
times nonafterburning thrust) at the expense of increased specific fuel consump- 
tion. The primary gas flow is then expanded through a nozzle to produce thrust. 
In some engines, the primary and secondary (bypass) gas flows are mixed ahead of 
the nozzle. These engines are called mixed flow exhaust engines. Other engines 
have separate nozzles for both the secondary and primary gas streams. Some 
turbofan engines have combustors located in the bypass duct and are, therefore, 
called ductburners or ductburning turbofans. 


METHOD OF ANALYSIS 


The method of analysis used in this study involves the thermodynamic modeling 
of each component (the compressor, burner, turbine, etc.) of a gas turbine 
engine. The analysis is applicable to a one or two spool turbojet or turbofan, 
with or without afterburning or ductburning. The analysis does not involve any 
cycle iterations. Properties of the combustion gases for use in the thermo- 
dynamic modeling are computed using polynomial curve fits of the specific heat, 
the enthalpy, and the entropy of the combustion gases (ref. 7). These fitted 
parameters are functions of the fuel -air ratio of the combustion gas and of the 
gas temperature. The technique is limited to temperatures below the dissocia- 
tion temperature (about 2200°K or approximately 4000°R). The properties in- 
cluded in the computer code are for JP-4, methane, or hydrogen; however, modi- 
fications to include any hydrogen, carbon, or hydrocarbon fuel are possible. 

The component efficiencies used in the thermodynamic modeling are scaled for 
off-design conditions from the input design-point value using empirical trends 
included in the code. 
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Principal Assumptions 


The principal assumptions in the present analysis are: 

1. Each compressor operates along a specified compressor operating line. 

2 . The bypass ratio is held constant at the design value. 

3. The flow through the turbine nozzles is sonic (choked). 

4. The combustor, afterburner and duct burner pressure drops are 
assumed constant at the design (input) value. 

5. The thermodynamic properties are calculated using polynomial curve 
fits (ref, 7) which include the effect of fuel addition, but do not 
include dissociation. 

6. The exhaust nozzles for the primary and secondary flow are assumed 
to be separate. 

A total enthal phy-entropy diagram for a two spool turbofan is presented in 
Figure 2. Each of the individual engine components are labeled and both the 
ideal process and the actual process are shown. Diagrams for other types of 
engines would be similar. 


GOVERNING EQUATIONS 


Gas Properties 

The gas properties have been approximated by the polynomial expressions of 
reference 7. These properties are functions of gas temperature, fuel to air 
ratio, and type of fuel. The computer code includes property descriptions for 
JP-4, hydrogen, or methane. 

Using this technique, the specific heat of the gas is expressed by a polynomial 
equation. Appropriate integrations of this expression yield the enthalpy 
change and the entropy function. The entropy function is defined as 


A(|) 



S dT 
T 


0 ) 
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The enthalpy change is defined as 



The importance of the entropy function can be seen from recalling that a 
general expression for entropy change between two state points is 



For an isentropic process this reduces to 


(4) 


When the appropriate integrations are accomplished, equations (1), (2), (3), 
and (4) produce a property description for the combustion gas. Using this 
property description, it is possible to compute the properties through the gas 
turbine from thermodynamic considerations. If the type of fuel is specified, 
the properties can be expressed in functional form. Using this form, the 
entropy function becomes 

A* = it f) 

the enthalpy function is 

I 

Ah = Ah (Tg, T^, f) 

and the specific heat is 

Cp = Cp (T^. Tp. f) 

where Tg is the initial temperature, T|^ is the final temperature and f is 

the fuel/air ratio for the considered mixture. Equations (5), (6), and (7) are 
also functions of the particular fuel being used (i.e. kerosene, hydrogen, 
methane). 


(5) 

( 6 ) 
(7) 
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Inlet 


The function of the inlet is to meter and condition the airflow to conform with 
the demands of the engine. During this process, there occurs a drop in total 
pressure of the flow because of losses due to friction as well as possible 
shock and separation losses. A measure of these losses is the total pressure 
recovery {the ratio of total pressure out of the inlet to total pressure into 
the inlet). The total pressure recovery will be used in this report to indicate 
the efficiency of the inlet. 

Standard functional relationships for inlet total pressure recovery have been 
developed (ref. 8), These usually relate the inlet total pressure recovery 
to the free stream flight Mach number. 

Two of these functional relationships are 

the AIA Ram Recovery 

Pl/Po = 1.0 > 0.1 (Mco - 1.0)^*^ M<»>1 

= 1.0 Moo<l 

and the Military Spec 5008B 

Pl/Po = 1.0 - 0.075 (Moo - 1.0)^ Moo>i 

= 1-0 Mco<l 

The AIA Ram Recovery and the Military Spec 5008B Recovery are included in the 
present program. It is also possible to input any other desired inlet total 
pressure recovery. 


(8) 

(9) 


Compressor 

Typical mathematical models of compressors involve detailed maps of compressor 
pressure ratio as a function of compressor corrected weight flow and corrected 
rotational speed with the compressor adiabatic efficiency as a parameter 
(refs. 9 and 10). These compressor models are usually empirically determined 
and are dependent on the particular compressor design. The empirical perform- 
ance of eight different compressor designs was examined in the article by 
Robbins and Dugan '(ref. 9). The purpose of their study was to obtain a general 
empirical description of compressor performance for a wide family of axial flow 
compressors. The present work uses the results of the study by Robbins and 
Dugan. In addition, the present -work utilizes the assumption that the compres- 
sor is operated along a line of .maximum adiabatic efficiency as the compressor 
weight flow and speed is varied. It is also assumed that the compressor refer- 
ence point described by Robbins and Dugan is equal to the compressor design 
point. This assumption eliminates an iteration and seems to cause no appre- 
ciable error. Using these assumptions and the work of Reference 9, it is 
possible to describe the performance of a family of axial flow compressors using 
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three charts which have been reproduced from Reference 9 as Figures 3, 4, and 5. 
Figure 3 presents the percent corrected speed as a function of design compressor 
pressure ratio and percent corrected airflow. The compressor corrected weight 
flow and design compressor pressure ratio are used to obtain the percent cor- 
rected speed using this figure. Figure 4 presents off -design compressor pressure 
ratio (in percent of design) as a function of design compressor pressure ratio 
and percent design corrected speed. The input to this figure is design compressor 
pressure ratio and percent design corrected speed which is obtained from the 
previous figure. The off-design compressor pressure ratio is obtained from this 
figure. Figure 5 presents off-design adiabatic efficiency as a function of design 
compressor pressure ratio and percent design corrected speed. The off-design 
adiabatic efficiency is obtained from this figure. Each of these charts is 
curve fit. If the cycle has multiple compressors, the above model is used on 
each separate compressor. 


After the pressure ratio and adiabatic efficiency for each compressor has been 
obtained, the work required to drive these components can be computed. The 
ideal compressor exit temperature can be found from 


A(J) (T 


out^. 




/P 


out' 


in 


(10) 


where Tg^^^_ is the 
The outer compressor 


ideal exit total temperature for the 
work is 





W. ^out- 
in\ 1 


- N-n 


outer compressor. 


( 11 ) 


where the compressor adiabatic efficiency is obtained from Figure 7. The inner 
compressor work is found in a similar manner. 


The total mass flow into the engine is the sum of the secondary (bypass) mass 


flow (Wgg(.) and the primary mass flow (Wpj-im) 


W = W + W • 

ov sec pnm 

(12) 

and the bypass ratio is defined as 


~ '^sec'^''*prim 

(13) 


The exit total temperature for the compressor can be found for the inner and 
outer compressor by using enthalpy obtained in equations (10) and (11), 
respectively, along with the definition of the incremental enthalpy in equation 
( 6 ). 
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Fan 


For bypass engines, the fan supercharges the primary flow. For the purpose of 
the paper, the inner portion of the fan flow which pressurizes the primary mass 
flow is considered part of the outer compressor. The outer annulus portion of 
the flow which enters the bypass duct, if the engine is a turbofan, will be 
denoted as secondary flow and any reference to fan performance will consist only 
of this outer portion of the fan. The fan and outer compressor are assumed to 
have the same pressure ratio per stage for both design and off-design operation. 
The fan off-design pressure ratio can therefore be determined by scaling its 
value from the outer compressor off-design pressure ratio by using the known 
values for the design fan and outer compressor pressure ratios. 


This scaling results in the equation 



/fan 


In (PRFD) 

^VTrTTFmj 

^in /comp 


(14) 


where p is the total pressure at the indicated points, and PRFD and PRCD are 
the design fan and compressor pressure ratios. 

Equation (14) can be used to approximate the fan pressure ratio at off-design 
conditions from the compressor off-design pressure ratio, which is obtained from 
the compressor model, and the ratio of design fan pressure ratio to design 
compressor pressure ratio. 

This technique eliminates the need to do cycle iterations to find the fan off- 
design pressure ratio. However, it is not as accurate as the techniques 
involving cycle iteration using separate fan maps and, in general, it over- 
predicts the off-design fan pressure ratio. The cycle iteration approach 
involves guessing values of the parameter of interest, computing through the 
cycle, checking the results against some required bounary condition, and then 
correcting the input parameter if the boundary conditions are not satisfied 
(refs. 1 to 5). 

The bypass ratio, the ratio of secondary airflow through the fan duct to primary 
airflow through the core gas turbine, is assumed to be constant at the design 
point value. 

The exit fan ideal total temperature can be found using equation (4) for the 
known fan pressure ratio. 

The adiabatic efficiency of the fan is assumed equal to the outer compressor 
adiabatic efficiency. The fan work can be obtained from an equation similar to 
equation (11). 


10 



ORIGINAL PAGE IS 
OP POOR QUALTPY 

Combustor/ Afterburner/Ductburner 


The computation for the combustor, afterburner, and ductburner involves a mass and 
energy balance across each of these components. If the mass balance and energy 
balance are combined, the following equation for the fuel -air ratio is obtained: 


W 


f = 


fin 


- h. 


W 


air • 
in 


product ~air 
(hf + Hg WF} - hproduct 


(15) 


The total enthalphy of the products (h , ) can be computed from the input 

products 

combustor exit temperature and the fuel -air ratio. The total enthalphy of the 
fuel is the sum of the sensible enthalphy (h^) and the chemical enthalphy 

(HVF) or lower heat value of the fuel. The program operation requires that the 
combustor design exit temperature be input. The combustor inlet temperature is 
reduced if the engine is to be operated at reduced power. 

The combustor pressure drop is assumed to be constant at the design (input) 
value. 

Typical values for the burner efficiency, On, are shown in figure 6 (refs. 11, 
12, and 13). The burner efficiency is shown as a function of P5 Ts/Vr, where 
Pj is the static pressure at the entrance, T5 is the static temperature 'at 

the burner entrance, and Vr is the combustor reference velocity, based on p^, 
T5, and the combustor maximum cross section area. The computer program contains 
a nondimensional version of this figure which approximates the lower limits of 
figure 6. 


Turbine 

The inner turbine usually sets the limit for the maximum combustor exit 
temperature. The inner turbine nozzles operate at a Mach number of one (choked) 
and thus limit the mass flow through the engine. For choked flow, it is assumed 
(ref. 14) that 

= constant (16) 

where W is the mass flow into the flow area A at a total temperature T 
and a total pressure p. In the computer program there are two modes 
of operation for the inner turbine compuxer model. 

The first mode is the design mode in which the turbine nozzle area required to 
supply airflow sufficient to produce an input thrust level at a given flight 
condition is determined. 




In this design mode of operation, the inner turbine nozzle flow area is computed 
for input component efficiencies, design pressure ratios and design combustor 
exit temperatures . The airflow through the engine is arbitrarily chosen for this 
mode as the airflow and flow areas scale directly with thrust and, therefore, the 
computed inner nozzle flow area can be scaled to the required input thrust level. 

The inner turbine nozzle area is then fixed for the second mode of operation - 
the off-design mode - and the combustor exit temperature is computed which satis- 
fies equation (16) for a known inner turbine nozzle flow area, mass flow, and 
total pressure. 


The isentropic turbine work is equal to the work required to drive the attached 
compressor plus the additional work required to overcome the losses of the 
turbine. These turbine losses are usually accounted for by using a turbine 
adiabatic efficiency, which is empirically determined. Thus, the isentropic 
work of the turbine is: 


isentropic 

turbine 

v/ork 


^comp out j n) 


comp 




(17) 


The isentropic total temperature out of the turbine can be computed from the 
ideal turbine enthalphy difference of equation (17) by use of equation (6). 

The turbine total pressure ratio can be computed from the turbine entry 
temperature, the isentropic turbine exit temperature and the fuel air ratio by 
use of equation (4). 


In order to complete the calculation, it is necessary to find the turbine 
adiabatic efficiency. The design point turbine adiabatic efficiency is input 
by the program user. The off-design turbine adiabatic efficiency can be 
computed from an approximation suggested by- Dugan in reference 9. The 
approximation can be expressed as: 


’^t " ’^t 

^design 



'-t. _ 1 


9 

' “ -1 
1 


[m! 

1 

in 

designj 


(18) 


Where q, . is the input design adiabatic turbine efficiency, N4. is the 
^design l 

turbine RPM and Ah is the actual total enthalphy difference across the 

turbine. The turbine RPM, N^, is equal to the associated compressor RPM. 

The total enthalphy difference across the turbine Ah can be found using: 


Ah 




(19) 
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Dugan stated that the accuracy of equation (18) for the off-design adiabatic 
efficiency was "quite good over much of the turbine performance map. Near 
limiting loadings, however, actual turbine efficiencies are much lower than 
those calculated." 

The outer turbine's performance is computed in a similar manner except the 
outer turbine flow passages are not sized or assumed choked. 


Nozzl e 


The nozzle calculations utilizes the assumption that the flow is fully expanded 
to ambient pressure isentropyically. A correction term, the nozzle velocity 
coefficient, is applied to the computed velocity to account for incomplete 
expansion as well as other losses. An energy balance across the nozzle for an 
ideal nozzle gives: 


^in " ^out “ *^out^ 



( 20 ) 


Using the isentropic fully expanded nozzle assumption equation (4) can be 
wri tten : 


nozzle 


R 

JM 



(21) 


Where Pout/Pin equal to the reciprocal of the engine nozzle total pressure 
ratio. The above expression for A(|) can be used to compute the ideal total 
temperature at the nozzle exit by using equation (5). The ideal nozzle total 
enthalphy can be computed using equation (6), and the ideal velocity can be 
computed from equation i[20). This velocity is then corrected to the actual 
jet velocity using an input velocity coefficient which can be a function of 
the engine nozzle pressure ratio (total pressure at nozzle entrance divided by 
ambient static pressure). 


Specific Thrust 

The specific thrust of an engine is defined as the net thrust divided by the 
total airflow of the engine,. The net thrust of the engine can be obtained 
from a momentum balance across the engine. Recall the nozzles are assumed to 
be fully expanded (Pexit ” ^ambient)- customary (refs. 8, 10, and 15) 

to reference the engine thrust to free stream conditions and make the appro- 
priate corrections for any external cowl drag or change in stream momentum 
ahead of the engine. Therefore, the engine thrust is 

Gross thrust = (Wg Vg + Wg Vg)/g 

Ram Drag = Wqv x Vg/g (22) 

Net Thrust = gross thrust - ram drag 

V 
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and the specific thrust using (22) is 


c-T net thrust 
“ W 

0 


''8 " ^6 ''6 - »n Vn 
“o 8 


(23) 


Specific Fuel Consumption 

The specific fuel consumption is defined as the total fuel flow in pounds per 
hour divided by the net thrust in pounds. This important parameter is a 
measure of the overall efficiency of the engine. The fuel flow per hour is 


fuel flow per hour = 





3600 


Using (22) and (24) with the definition of SFC 


SFC 


I "8 “8 " «6 ^6 - tJo M 


(24) 


(25) 


Overall Engine Efficiency 

The overall efficiency of the engine is another important parameter of the 
engine cycle. This efficiency is described quite completely by Shepard 
(ref. 8). The overall efficiency is defined as the useful power (thrust x 
velocity) divided by the net energy input per unit time. 

(net thrust) (Vo) 

n = 

ov rate of energy in 


The rate of energy into the engine is equal to the fuel flow times the energy 
stored in the fuel. Using the definition of SFC and considering both the 
chemical, and kinetic energy of the fuel, the rate of energy is 


rate of energy in 


SFC(net thrust) (HVF)«3 + 
3600 


(27) 


3600(Vjj) 



^71 

SFC 

^(HVF)J + 2g 1 


( 28 ) 
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In many texts, the kinetic energy of the fuel is neglected and thus (28) 
becomes 


3600 (V^) 

%v * (SFC)(HVF)(0) 


(29) 


The overall efficiency, can be divided into two components; the propulsive 
efficiency, and the thermal efficiency The propulsive efficiency 

is the useful power (thrust times velocity) divided by the rate of kinetic 
energy production by the exhaust 


(net thrust) (V q) 


P>-0 ^ (v/ . v/) . ^ . yi ) 


(30) 


The fuel added to the engine has been neglected in (30). The propulsive 
efficiency is a measure of the losses due to difference in exhaust jet velocity 
and the free-stream velocity. The thermal efficiency is the energy rate out 
of the engine which is available for useful work divided by the rate of energy 
into the engine. This is a measure of the efficiency of the cycle. 

Therefore, the thermal efficiency is 


W. 


\h ‘ lSFC)ln . |t^th . r u s ii t v//2g] 

and note that by comparing (28), (30), and (31) 


(31) 


^ov “ ^th ^pro 


(32) 


ACCURACY OF THE PROGRAM 


The accuracy of the program can be evaluated by comparing the program output 
for a known engine cycle with the comparable specification (ref. 16). of the 
known engine. This comparison requires the knowledge of the component effi- 
ciencies for the engines which are used in the comparison. Unfortunately, the 
component efficiencies are not available in the literature for particular 
engines, however, general trends and typical values are well known (refs. 9, 
10, 15, and 17) and could be used for the comparison. Another method of 
obtaining values for the component efficiencies for a particular engine is to 
match the sea level static performance of the engine by varying the input 
component efficiencies within the bounds of known technology. The approach, 
using this technique, is the fix all known engine cycle parameters and to vary 



the component efficiencies until the sea level static performance from the 
cycle analysis is equal to the actual engine specifications. The results 
using this technique are not unique but experience shows that values for the 
engine component efficiencies do not vary more than a few percent from design 
to design if the engines were designed at the same technology level. In the 
following comparisons, the sea level static 100 percent RPM, thrust and SFC 
point was matched approximately to the engine specifications (ref. 16) by 
varying the component efficiencies. These efficiencies were examined to be 
sure that they were reasonable and consistent. The code then was used to 
generate all off-design (part throttle) and design points other than the sea 
level static, 100 percent RPM point. During these calculations, the input 
design point efficiencies were used, with off-design variations to these 
parameters computed from curve fits of empirical data contained in the code. 
Figures 7, 8, and 9 present comparisons for the JT8D, the J-79, and the CF6 
engines, respectively. The engine specifications from Reference 16 are repre- 
sented by the solid lines. Average differences of 5 percent are typical with 
maximum error of less than 15 percent. The maximum error usually occurs 
during the after-burning or duct-burning and is probably caused by the assump- 
tion of constant combustion efficiency and pressure drop which was made in the 
code. 


RESULTS AND CONCLUSIONS 


A method for computing the design and off-design performance of a gas turbine 
engine is presented and the results are compared with engine specifications 
for several typical engines. The results indicate average errors in the 
performance calculation of less than 5 percent with maximum errors of about 15 
percent. The maximum errors usually occur when the after-burner is in 
operation. 

A computer code is described which is capable of producing performance 
calculations for a gas turbine engine at one Mach altitude and power setting 
in less than 0.1 second on the IBM 360/67 as compared to slightly less than 10 
seconds for GENENG* 

The increase in computational speed was accomplished by not using detailed 
compressor and turbine maps and not using any iterations between components in 
the calculations. The compressor model us.ed in the present report is based on 
empirical data and required the assumption that the compressor Is operated 
along a line of maximum compressor adiabatic efficiency. 
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Figure 1.- Engine schematic for a two spool turbofan. 
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APPENDIX A 


Program Input 


Input to the computer program is achieved via a namelist called LEWIS . 
The various input variables are described in this appendix. The inputs 
followed by ** are required for operation of the program. The other 
Inputs have default values incorporated in the code. 
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APPENDIX B 


Program Output 

1) DETAILED OUTPUT NOMENCLATURE 

2) SAMPLE DETAILED OUTPUT (PWJT8D-9) 

3) SAMPLE SUMMARY OUTPUT (PWJT8D-9) 



DETAILED OUTPUT NOMENCLATURE 




pwrc 

PPPCFMT OFc:TOM CnPPFCTFD AIRFI.OW 


^PK'THF 

PPPCFMT OFFTOM CORRKFTFD FAM RPFFO 

4 

PMTHC 

PFPCF^!T OFFTON CORRECTPO HIGH PPFSS COMPRESSOR 

SPFFO 

N TPPUST 

MFT Twpu«;t (LRF) 


OP THRUST 

GROSS THPIIST (LRR) 


RAW nPAG 

RAM DRAG (I.PF) 


orr 

SDFriFTC FUF| COMSt'MPTTON (LRm/hR/LRF) 


PROP FFF 

PRODIILSI'/R FFFTCIEfjCY 


OVl FFF 

nVPRAU. frftctfucy 

* 

WTOT 

total MARSFI OW PFR PNGTMF (LRM/SFC) 


WCOPR 

rORRFrTFO wASSFLOW PFP ENGir^F (LRM/SFC) 


SPFF THR 

SRFCTFTC TMPDST ( LBF/ (L8H/SFC) ) 


FUFL flow 

FUFl FI OW PFP FN'GIMF (! BM/HR) 


WCV’APRT 

PATTO OF ACCOUNT AHLF 'COOLIM^ MAqSR.OW TO TOTAL 

PRIMARY MASSFLOW 

WRLFVPRI 

P&TTO OF RLFFO vASSFLO-i'-TO- 'TOTAL RPTMARY MASSFLOW 

POWFR 

PO>''FR FXTPACTFD from FMGTNE tPTU/SFC) A POSITIVE NUMBER 

DPS RR 

OFSTGM BYPASS RATIO (OtlCT TO PRIMARY FtOW) 

• 

DE . OPR 

OESTGN overall COMPRFSSOP -'PRFRSURF RATIO 


DES FPR 

OFSTGN FAM PRESSURE PATIO 


dfr ttt 

design COMF1 'STOP EXIT TEMPERATUPF '(OF-G R) 


ACT PR 
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SI 
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3 k 


ACT OPR 


actual OV'TRALL CQMPQFSSOR ppp^supf ratio 












SAMPLE DETAILED OUTPUT 
PW JT8D-9 




I-TOT = t.'rOP9 = THP = <MS,nOQ FUEL FLOW= Fftg?, 


„C'>‘Ai^'PT = .ripp vo(_r.[>f^T = n.uon poi-F'^ = 0,0 

rtc; up = i.r-io pp<; o:^!? s: PFS FPR = 1,^10 HF^: TIT =PIPO.0n 
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F19.7 

Sir, 7 

“ FArr _ 

pp = 1,010 fTC = ,8Pon P(\iTH = l(>n, 000 
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ni^V T(iOv>0C/^,ri FKiftilMF rYCl.F pp 
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^ THD|Jc;t = 1437.0 r.^ TH‘">nsTs 33SS,‘=, tjAM pR/iG 

OFC -» ,7^47 PPOP FFF = .7^»^Q QVL PpF 

WTOT s 7F.PP wrriPP =r 30?, R<? $PFr TNQ 

^■' C>^'flP»^?T - .0?0 = O.nOQ PnwFR 

ppo nc = l.nin nFF opo poP 

ACT flP = 1.0’=*0 ACT QPP = I'^.OPP ACT FPP' 


GPAM 


s: 06,00 

PWTHF = 

Q7.44 PKJTHC 

s 07,40 

= 1P17,7 
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• - 

= ,?6H3 
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FUFL FLOWS 
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0.0 


• 



= 1 .910 nF«^ ttt sPiRo.nn 


l.Bh? ACT TIT =1617,60 












VACH = .p*0 Al.TrTtinp- =/.M)OO.P Pvrc = P4.00 PMTHF « 9fy,JA PNTHC' = 96. n 


M THPn«:T = 1043.*^ riP PAM fWAR s: 1H77.7 

«;fC ■= ,7<}c,? PPOP c-FF 5= .7^40 0V|. FFF = ..P646 



74.^0 '-'CODO = pqfi,^cjq <iPFC THR = l4,Sft9 FUFL FLOW* ft6?. 
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SAMPLE SUMMARY OUTPUT 
PW JT8D-9 




"tHniiQT' 


T IWFFT 



















.34 

35000, 0 

.96.0 

60?4.H 

3156,3 

?167,7 



,64 

35000.0 

96.0 

5416.0 

3091.9 

1766.1 


, 


35000.0 

94.0 

4T6?,5 

3ti?7.5 

1403,1 


• 


:^'=^r(On,n 1076.4 


.«4 36ono,fl <?n.n 316?. 6 ?rvh,7 


,H'+ 4'^ono.oi 00.0 40U^,4 IVQn.n 1604,0 


46000.0 «}R.n 3730,1 1966.1 133?. 0 
i 4SPO0.0 96,0 3364.0 1916.1 1064,1 
4 440^0,0 94,0 pq^n.? 1476,? 36?. 0 
4' 450 00.0 9?7n ?5?4,3 661 . \ 
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